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Abstract 

In this paper we are concerned with the problem of a crack lying at the interface between dissimilar 
materials with microstructure undergoing antiplane deformations. The micropolar behaviour of the 
materials is described by the theory of couple-stress elasticity developed by Koiter (1964). This 
constitutive model includes the characteristic lengths in bending and torsion and thus it is able to 
account for the underlying microstructure of the two materials. We perform an asymptotic analysis 
to investigate the behaviour of the solution near the crack tip. It turns out that the stress singularity 
at the crack tip is strongly influenced by the microstructural parameters and it may or may not show 
oscillatory behaviour depending on the ratio between the characteristic lengths. 
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Nomenclat ure 



G± 


shear modulus (+/— stands for upper/lower half-plane) 




Poisson's ratio 


, + 


material characteristic length in bending 




material characteristic length in torsion 


h 


material characteristic length (Koiter's notation) 


n± 


ratio between the characteristic lengths in bending and torsion (Koiter's notation, — 1 < t)± < 1 


u 


displacement field 


w 


out-of-plane displacement 


V 


rotation vector 


e 


strain tensor 


X 


curvature tensor 


t 


nonsymmetric stress tensor 


a 


symmetric part of the stress tensor 


T 


skew-symmetric part of the stress tensor 




couple-stress tensor 


P 


reduced tractions vector 


q 


reduced couple-stress tractions vector 


A 


order of stress singularity 
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1 Introduction 



Nowadays, bimaterials are efficiently and widely used in many advanced engineering applications, such 
as layered composite structures, electronic packaging and thin film coatings. For the prediction of failure 
of these structures and the assessment of acceptable stress level under the condition experienced during 
service, it becomes essential to estimate the magnitude and distribution of the interfacial stress and 
strain fields along the interface and mainly near the tip of interface cracks, which may arise and extend 
under general loading conditions. In particular, antiplane shear loading condition may frequently occur 
in the life span of composite structures, both alone or accompanied by plane deformation. 

Within the classical LEFM theory, the crack tip fields for an interface crack under antiplane strain 
are similar to the Mode III crack tip fields in a homogeneous medium (Willis, 1971; Piccolroaz et al., 
2009). In both cases, indeed, the shear stresses on the crack plane are the same in the upper and 
lower bodies, whereas the out-of-plane displacement is zero on the uncracked region of the crack plane. 
Thus, it is possible to combine the lower and the upper bodies to obtain equilibrium, without changing 
displacements or stresses in the two halves. Stresses exhibit Mode III symmetry, but displacements do 
not and thus crack sliding profiles are not symmetric. 

Due to the lack of a length scale, the classical theory of elasticity is not able to characterize the 
constitutive behaviour of brittle materials at the micron scale. This lack is expected to be particularly 
significant for the analysis of the stress and deformation fields very near the crack tip. For a proper 
investigations of the crack tip fields at the micron scale it becomes necessary to adopt enhanced con- 
stitutive models, which account for the presence of microstructure. A way of doing that consists in the 
inclusion of one or more characteristic lengths, typically of the same order of the compositional grain 
size, generally few microns, for many advanced materials. The indeterminate theory of couple-stress 
elasticity (CSE) developed by Koiter (1964) involves the material characteristic lengths in bending and 
torsion. It is sufficiently accurate to simulate the behaviour of materials at the micron scale as well as 
the size effects occurring at distances to the crack tip comparable to characteristic lengths, but it is also 
simple enough to allow the achievement of closed- form solutions. 

Although the presence of the microstructure is expected to modify the interface crack tip field with 
respect to the classical solution of the LEFM, no analytical investigations have been so far performed 
about the problem of an antiplane crack along the interface between micropolar and classical elastic 
materials (the only related work regards a crack terminating perpendicular to the interface; Mishuris, 
1985). Most of related studies available in literature instead concern the problem of an interface crack 
under plane deformations, e.g., Itou (1991) examined the effect of couple-stresses on the strain energy 
release rate for an interface crack loaded by an internal pressure, neglecting somehow the oscillatory 
behaviour of the crack tip fields. 

In order to provide an experimental basis for studying the interfacial behaviour of a bimaterial 
specimen under shear loading, Kang et al. (2002) applied a method which combines moire interferometry 
with phase shift and image processing to measuring the interfacial displacement and strain fields within 
the interfacial region. Their experimental results show that there is a boundary layer characteristic with 
a peak value of shear strain and high gradient of rotation angle in the interfacial region. Their study also 
shows that similar results can be analytically predicted by means of the couple-stress theory, considering 
the additional freedom of the rotation angle effect. 

Hutapea et al. (2003) investigated the micro-stress generated along a fiber/matrix interface under 
generalized plane deformation, which are expected to dominate the failure initiation process in composite 
laminate. They showed that the micropolar theory is able to capture the interface micro-stress accurately. 

A small number of interface crack problems have been investigated by using the strain gradient theory 
of plasticity (Hao and Liu, 1999; Chen and Wang, 2002). In particular, Hao and Liu (1999) analyzed 
the crack propagation in bimaterial systems showing that high stress triaxiality always occurs on the 
softer material, which may promote ductile damage and facilitate crack growth. Chen and Wang (2002) 
explored the interface crack tip fields at micron scales under plane strain conditions. Their numerical 
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investigations show that the singularity of stresses in the strain gradient theory slightly exceeds or 
equals to the square-root singularity independently of the material hardening exponents. Askes and 
Gitman (2009) showed numerically that in gradient elasticity no singular behaviour is found for a crack 
terminating perpendicular to the interface. 

The analysis of singular stress concentration in homogeneous micropolar elastic solids shows that 
several pathological predictions of classical elasticity in singular stress concentration problems are altered, 
mitigated, or possibly eliminated when couple-stresses are taken into account (Nazarov and Semenov, 
1980). In particular, the problem of a Mode III crack in a homogeneous materials modelled by the 
couple-stress elastic theory was first analyzed by Zhang et al. (1998) and later by Geogiadis (2003) by 
considering a single characteristic length. The results obtained therein indicate that the skew-symmetric 
stress components have r -3 / 2 singularity near the crack tip, where r is the distance to the crack tip. 
Although this singularity is much stronger than the conventional square-root singularity, it does not 
violate the boundness of strain energy surrounding the crack tip and leads to a finite energy release 
rate. Their asymptotic analysis also provides a negative out-of-plane displacement ahead of the crack 
tip. This unphysical result is due to the exclusion of the lowest order terms for the displacement and 
symmetric stress components, which do not contribute to the energy release rate. 

The effects of both characteristic lengths in bending and torsion and a complete investigation of the 
crack tip fields under Mode III loading condition in homogeneous CSE materials have been properly 
addressed in a recent work by Radi (2008). The roles of both characteristic lengths are therein examined 
in detail and their influence on the crack tip is analytically explored by using Fourier transform and the 
Wiener-Hopf method. The asymptotic and full-field analyses show that the symmetric stress is finite 
at the crack tip, whereas the skew-symmetric stress is negative and strongly singular. Ahead of the 
crack tip within a zone smaller than the characteristic length in torsion, both the total shear stress and 
reduced tractions occur with the opposite sign with respect to the classical LEFM solution, as predicted 
by the asymptotic analysis. However, the zone of dominance of the asymptotic fields has limited physical 
relevance and becomes vanishing small for a characteristic length in torsion of zero. In this limit, the full- 
field solution recovers the classical Km field with square-root stress singularity. Outside this zone, the 
total shear stress exhibits a positive maximum, thus providing more realistic predictions on the tractions 
level ahead of the crack tip than the singular LEFM solution. A sharp crack profile is also observed. 
It may denote that the crack becomes stiffer, thus revealing that the presence of microstructures may 
shield the crack tip from fracture. 

In the present work, the effects of strain rotation gradients on a stationary antiplane crack along the 
interface between two different couple-stress elastic materials are analytically investigated by performing 
an asymptotic analysis of the crack-tip fields. The special problem of a crack along the interface between 
a couple-stress elastic solid and a classical elastic medium is also addressed in Sec. [3l The results of the 
present asymptotic analysis are expected to hold in a small zone near to the crack tip, whose extent 
may vary with the size of the characteristic lengths, and provide valuable information for performing a 
full- field analysis of the interface crack problem, e.g., by using the Wiener-Hopf method, which will be 
the subject of further investigations. 

2 Crack at the interface between couple-stress elastic materials 

We consider a bimaterial plane made of two dissimilar materials, joined along a perfect interface. The 
two materials are assumed to have an underlying microstructure, described by the material characteristic 
lengths in bending and in torsion, denoted by if and if , respectively. The elastic moduli are denoted by 
G± (shear modulus) and v± (Poisson's ratio). A semi- infinite plane crack is placed along the interface, 
and a Cartesian reference system is assumed centred at the crack tip, see Fig. [T] 

The fundamentals of the Couple Stress (CS) elasticity theory (Koiter, 1964) can be found in several 
text books and research papers (see for example Nowacki, 1985; Asaro and Lubarda, 2006: Radi, 2008). It 
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Figure 1: A crack at the interface between dissimilar couple-stress materials. 



is recalled here that the main characteristic of this theory is that the rotation vector ip is not independent 
of the displacement vector it, but it is subject to the condition 



1 



curl tt. 



(1) 



Consequently, all the kinematical quantities can be derived from the displacement field. In particular, for 
antiplane shear deformations, the following kinematical relations between the out-of-plane displacement 
w, rotation vector y>, strain tensor e and curvature tensor \ are derived 



1 dw 
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1 dw 
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-X22 = r 
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An infinitesimal surface element transmits a force and a couple vector, which give rise to a nonsym- 
metric stress tensor, t, and a couple-stress tensor, [i. The nonsymmetric stress tensor can be decomposed 
into a symmetric part er and a skew-symmetric part r, such that t = cr + r. 

The isotropic constitutive equations are given by 



2i/G(tre) 
1 - 2i/ 



/ + 2Ge, /x = 2G/ 2 (x T + TO), 



(3) 



where I is the identity tensor, I and the CS parameters introduced by Koiter (1964), with —1 < ?/ < 1, 
the superscript T denotes transposition. The material parameters I and 77 characterize the microstructure 
of the material and can be expressed in terms of the material characteristic lengths in bending and in 
torsion as follows: 

h = 1/V2, i t = 1^1+^. (4) 

For antiplane shear deformations, the nonzero stress and couple-stress components are 



dw 



dw 
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dx\dx2 



fJ-12 



dx\ 



ni 2 I d2w 



dx 2 ' 
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dx\ 
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d 2 w d 2 w s 



dx\ 



dx\ 
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In the absence of body forces and body couples, the equations of equilibrium read 



divt J 



0, div/x T + (rei) xejl (re 2 ) x e 2 + (re 3 ) x e 3 = 0, 



(7) 



where {ei, e2, 63} is an orthonormal basis. For antiplane shear deformations, the nonzero skew-symmetric 
stress components, derived from © and (O2, are 



T13 



Gl 2 dw 
~ ~dx~i 



T23 



Gl 2 dw 
~2~ dx~2~' 



(8) 



where A stands for the laplacian operator. 

In the CS theory, due to the internal constraint ([1]) between rotations and displacements, the Neu- 
mann boundary conditions are prescribed in terms of the so called reduced force tractions vector p and 
couple tractions vector q defined as 

1 



P 



t T n 



jV/Xnn x n, q 



T 

fi n 



(9) 



respectively, where n denotes the outward unit normal and fj, nn = n ■ \xn. Additionally, if the external 
surface of the body is not smooth but piecewise smooth, the boundary conditions include the following 
equation along each edge 

Q = o(/4n -Mnn)> ( 10 ) 



where Q is a prescribed line load tangential to the edge (Koiter, 1964), and superscripts + and — stand 
for the values on the surface at each side of the edge. It then appears that the condition (fT0| becomes 
essential in the case of bodies with non-regular boundaries, such as cusps, wedges and cracks. 

A substitution of §5§ and © into ©i gives the following governing equation for the displacements 



in the two half-planes: 



l 2 

2 



0, 



(ii) 



where A 2 denotes the bilaplacian operator. 

We assume that the crack faces are traction-free, so that the following boundary conditions apply 
for x 2 = and x\ < 0: 
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dx\ 
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(12) 



Assuming also that no tangential line load Q is applied along the crack edge, we enforce that 



M22 - -"22 = -G+Z+(l + ?7+) 



i:) 2 



8x18x2 



+ G-l 2 _{\+T 1 -) 



if 



(x u x 2 ) = (0- ,0+) 



dx\dx2 



0. (13) 



(zi,a:2) = (0-,0-) 



The formulation is completed by the transmission conditions for ideal interface, which imply con- 
tinuity of the displacements, rotations, reduced stress and couple-stress components for x 2 = and 
Xi > 0: 



M = 0, 



= 0, 



G 



( dw I 2 d 
\ dx 2 2 dx2 



(2 + V ) 



dx\ 



dx\ 



= 0, 



GV 



dx 2 



dx\ 



where the notation \f\ stands for the jump of the function / across the interface: J/J] 



= 0, 



(14) 
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2.1 Asymptotic analysis and singularity of stresses 

Assuming a polar reference system centered at the crack tip, we search for the main asymptotic term of 
the solution asr^O in the standard form as follows 

w ± (r,9) = r x F±(9,\). (15) 

We are interested in finding the leading term of the asymptotic solution corresponding to finite elastic 
energy. This requires that A > 3/2 (see Radi, 2008). It is noted that the values A = and A = 1 are 
also admissible, as long as the respective terms in (fl"5|) correspond to a rigid body motion (constant 
displacement) and a uniform deformation (linear displacement), respectively. Moreover, the expression 
(|15p can be used to find more terms in the asymptotic solution in the form ^ i r Xi F± (9, Aj), provided 
that |A, — Xj\ < 2, V i ^ j. If more terms are required with exponents differing by 2 or more than 2, 
then a two-terms asymptotic procedure should be used instead, as explained in Sec. [3] 

Keeping into account only the leading term as r — > 0, the governing equation ([lip yields the following 
ODE for the unknown functions F± 

F±" + 2(A 2 - 2A + 2)F'l + A 2 (A - 2) 2 F± = 0. (16) 

We first investigate the simplest cases A = and A = 1, for which eq. (|T6|) admits the solutions 

F±(6, 0) = Bf + Bf9 + Bf sin 2(9 + Bf cos 261, (17) 

Ft (0, 1) = ( B t + B 2°) sin 9 + ( B 3 + B A e ) cos 6\ ( 18 ) 

respectively. Taking into account all boundary and transmission conditions, one can conclude that for 
A = and A = 1, eqs. (TIT)) and (TT8"|) take, as expected, the forms 

F±(0,O)=/3 O , f±(M) =Ai sine + /3 12 cos 9, (19) 

respectively. 

One can also use the representation (|15|) to find the solution for the case A = 2, since the term 
corresponding to A = is not involved in the analysis (as it vanishes after differentiation). Thus, for 
A = 2, eq. (TlB]) admits the solution 

F±(0,2) = Bf + Bf9 + Bf sm29 + Bf cos28, (20) 

which, taking into account all boundary and transmission conditions, reduces to 

F± (9,2) = (3 21 (cos 2 9 + ? ]± sin 2 9) + f3 22 sin 29. (21) 

For all other cases, eq. (fT6|) admits the following solution: 

F± (9) = Bf sin(A6») + Bf cos(A6>) + Bf sin[(A - 2)9] + Bf cos[(A -2)9}. (22) 

By imposing the boundary and transmission conditions, we obtain a 8x8 homogeneous algebraic system, 
whose characteristic equation is 

sin 2 (7rA) [cos(2ttA) + k) = 0, (23) 
where k = C/D > for any r/± > — 1 and 

C = G\l\{b - 2ry_ + V 2 _ )(3 - r/+) 2 (l + r /+ ) 2 + Gilt (5 - 2 V+ + ^ )(3 - ?7_) 2 (1 + r?-) 2 
+2G + l 2 + G-l 2 _ (3 + 77+ + v- - - »7+)(l + r?+)(3 - + r/_), 

(24) 

D = (3 - r?+)(l + r?+)(3 - r?_)(l + r?_) {G^|(3 - r? + )(l + V+ ) + G 2 _lt(3 - r?_)(l + r?_) 
+2G+Z 2 r G_Z 2 (5 - r?+ - r/_ + r?+?7_)} . 



The first term in eq. (|2lt)) leads to the conclusion that A = 3 (the cases A = 0,1,2 have been 
investigated above), while the second term may exhibit singular behaviour depending on the value of 
the parameter n. If k > 1, then the solution of eq. (|23p is complex and the singularity shows oscillatory 
behaviour in the vicinity of the crack tip. Otherwise, the solution is real and there are no oscillations. 
More precisely, since k is strictly positive, the following three cases may occur: 

(i) < k < 1: the first admissible value of the exponent is in the interval 3/2 < A < 7/4 (simple 
root). 

(ii) k = 1: the first admissible value of the exponent is A = 3/2 (double root). 

(iii) k > 1: the first admissible value of the exponent is A = 3/2 ± 17 (simple root), where 



7 = — log( K +^ 2 -l)- (25) 



In the case of a homogeneous material, G+ = G-, 1+ = 1-, i]+ = V-> the ratio k is equal to 1 and 
thus the first admissible value for the exponent is 3/2 (this case has been analysed in Radi, 2008). Some 
other special cases are investigated in the next section. 



2.2 Some particular and special cases 

To decrease number of parameters, let us first consider the case where 77+ = ?/_ = 77. Then the ratio k 
reduces to 

= {G\l\ + G 2 _l 4 _)(5 - 2?? + rj 2 ) + 2G+l 2 + G-l 2 _ (3 - 77) (1 + rj) = (a 2 + b 2 )c + 2ab 
K ~ {G\l% + G 2 _li){3 - 77XI + 77) + 2G+l 2 + G-l 2 _{b - 277 + r? 2 ) ~ a 2 + b 2 + 2abc ' ( ' 

where we use the notations a = G+l\, b = G-l 2 _ and 

5 - 277 + V 2 

C{V)= (3 -„)(! + „)■ (2?) 

Since c > 1 for any admissible value of 77 (—1 < 77 < 1), it is easy to show that, for dissimilar 
materials, k is always greater than 1 and equal to 1 if and only if c = 1 (or equivalently 77 = 1), see Fig. 

m 

In the limiting case 77 = 1 the asymptotic solution in the vicinity of the crack tip is given by 



I 2 ( ( 9 39 
w ± (r, 9) = p + rtfi sin 9 + fa cos 9) + r 3 / 2 -g- j fa ( 3 sin - - sin — 

( 6 39 \ 1 

cos --cosy J \ + r 2 (fa + fa sin29) + 0(r 5/2 ), r -> 0, 



(28) 



where the constants fa are amplitude factors depending on far-field loading and specimen geometry. It 
is noted that logarithmic terms are excluded, since geometric and algebraic multiplicities of the double 
root coincide (Rossle and Sandig, 1996). 

Correspondingly, the asymptotics of symmetric stress, couple-stress and skew-symmetric stress as 
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Figure 2: Plot of the ratio k given by l26l as a function of r\ for a/6 = 1, 3, 6, 10. 



are 

af 3 = f3 2 G ± - 2r x '% sin 2 t (^3ft sin °- + ft cos + 2rG± (ft cos 9 + ft sin 0) + 0(r 3 / 2 ), 

of 3 = ftG± + r 1 / 2 / 2 , sin - (3ft sin 6» + ft (3 + cos 0)) + 2rG± (ft sin + ft cos 0) + 0(r 3/2 ), 

^ = -r- 1/2 l 2 + Z 2 sin - [3ft (sin + sin 20) + ft (2 + cos0 + cos 2(9)] + 4ftG±Z± + 0(r 1/2 ), 

M21 = /-*12 = r" 1/2 Z^ Z 2 ^ I ft cos I (sin 26> - sin 0) + ft cos |(2 - cos0 + cos20) j + 0(r 1/2 ), 

rt 3 = y s/2 l 2 + l 2 - (3ft sin I +ftcos|) + 0^), 
1 Q/o 9 9 . 3$ 3$. 



(29) 



_± 



To 



= --r-VH 2 + l 2 _ (3ft cos — - ft sin — ) + O^ 1 / 2 ). 



Applying now the condition (j 13|) . we obtain 

A1 +(0 = 7 r)- Al22 (^-7r) = ^±^+4ft(G_/ 2 -G+Z^ + O^ 1 / 2 ), r -> 0. (30) 

v r 

The limit, as r — > 0, of ((50)) should equal the tangential line load Q applied to the crack edge. It then 
appears that ft is always zero and ft does not vanish only if Q is different from zero. It is also noted 
that the constant ft plays the role of stress intensity factor. 

The asymptotics for the crack opening, [w]] = w + (9 = n) — w~(6 = —it), and for the skew-symmetric 
stress ahead of the crack tip, t^(0 = 0), as r — > 0, are given by 

M = 4ft G+ ^ + + /^- r 3 / 2 + 0(r 5 / 2 ), t±(0 = 0) = r" 3 / 2 + O^ 1 / 2 ). (31) 

Therefore, in the vicinity of the crack tip, the skew-symmetric stress at 9 — displays a sign opposite 
to that of the crack opening, in contrast to the classical result of LEFM. A similar effect has been found 
by Radi (2008) in the case of a crack in an homogeneous CS material. 
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In the opposite case (77+ = = 77 < 1), the solution always exhibits oscillatory behaviour near the 
crack tip. Moreover, this behaviour is quite different from that we encounter in the case of classical 
materials. In the classic case, the region near the crack tip where the oscillatory behaviour appears is 
very small, while in the considered case this zone can be quite pronounced or its size can even tend to 
infinity if 77 — > — 1. Taking into account that the asymptotic analysis given here is valid only in a small 
neighbourhood of the crack tip where the micropolar theory controls the behaviour of the solution, such 
situation has limited physical meaning. 

Let us now assume that the parameters 77+ and rj- are different and one of them, say tends to 
the limiting value 77 _ — > —1, while the other is separated from — 1, so that 77+ > — 1 + e, where e is 
a small positive parameter. In this case, one can easily check that k — > 00 for any fixed value of 77+. 
This again corresponds to the case where the exponent is a complex number with the imaginary part 
becoming infinite, 7 — > 00, and therefore it has no physical relevance. 

It is noted that the case r?_ = —1, 77+ 7^ —1 cannot be recovered from the limiting case discussed 
above, and thus it will be discussed separately in Sec. 13.31 

Another case of interest is when one of the multipliers involved in the parameter k, say G-P_, 
diminishes, G-lt_ — > 0, then 

5 — 2ti_ + 7? 2 

^ rara =cy - L (32) 

Once again, this solution has physical relevance if and only if r/ = 1. The assumption G—l— — > takes 

place if one assumes that /_ — > 0, or, in other words, when the material in the lower half-plane reduces 
to a classical elastic material. This suggests considering the problem of a crack at the interface between 
a micropolar material (occupying the upper half-plane) and a classical one (occupying the lower half- 
plane). However, it is not possible to recover from ([28]) the solution for a classical elastic material, due 
to the singular perturbation character of the equation (fTTj) as l± — > 0. For this reason, the problem of 
a crack lying at the interface between couple-stress elastic and classical elastic materials is addressed 
separately in the next section. 

2.3 Energy release rate 

In this section, the energy release rate is evaluated for the asymptotic representation (f2"5|) . valid in the 
case ?7_|_ = 77- = 1, by means of the J-integral argument. The conservation law for couple-stress elasticity 
(see Lubarda and Markenscoff, 2000) implies that 

J= f ( Wm - t T n ■ |^ - n T n • ds = 0, (33) 

Jr V 9xi 9xi J 

for every closed contour T, provided that there is no singularity within T. In (|33p . n is an outward unit 
normal on T and W denotes the strain-energy density 

W = Ge-e + Gl 2 ( X -X + VX-X T )- (34) 
We define r± = Tf U Tf U T± U (see Fig. [3J, so that 

^(J 2 ± TJ±-J 1 ± ±4) = 0, (35) 
± 

with evident meaning of the symbols. Then we evaluate each term according to the representation (1281) . 
thus obtaining 

, 2/3 2 /3 4 ^ 2 , 1/2, , s 

Ji = Ji +Ji=- %± - + P + 0{R\'\ (36) 

V-fll 
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-+P + 0(R 1 2 /2 ), 



15 



271 



±- 



V^2 



0(Rl /2 ), 



+ 0(Rl /2 ), 



in which 



(37) 
(38) 
(39) 

(40) 



Note that the integrals along the interface cancel out from whereas the integrals along the crack 
faces give 

J--J+=F(R 2 )-F(R 1 ), (41) 



where 



F(R) = 



R 



0(R 



(42) 



It is possible now to introduce a path- independent parameter as follows: 



J* = Jim (h + F(R 2 )) = lim (Ji + F(Ri)) = P. 



(43) 



In consideration of the condition (|13[) , (3^ vanishes and the path- independent parameter J„ gives the 
energy release rate (derived through direct energy balance considerations by Atkinson and Leppington, 
1977 and Eshelby, 1980) 



J* = yKl+l- 



ft 



(44) 
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In the case of homogeneous body (G+ = G_ — G, 1+ — I- = I), our formula coincides with earlier 
results by Radi (2008): 

J* = ^l (45) 

It is noted here that particular attention should be paid when using the J-integral in CS materials. 
The contribution of the integrals along the crack faces is not zero, since the force tractions t T n and 
couple tractions fi T n are not vanishing as in the classical case. 

3 Crack at the interface between couple-stress elastic and clas- 
sical elastic materials 

In this section we consider the problem of a crack lying at the interface between a micropolar material 
(occupying the upper half-plane) and a classical elastic one (occupying the lower half-plane). Then the 
governing equations are 

p 

Aw + - y A2w+ = °' Aw ~ = °- ( 46 ) 
The traction-free boundary conditions along the crack faces are: 

p+ = and q+ = for 6 = ir, ovT 3 = for 6 = -tt. (47) 

Along the ideal interface the continuity of displacement and force tractions needs to be enforced, so 
that we have the following transmission conditions: 

w+ = w~ , pj = CT^ for = 0. (48) 

However, since the orders of the two governing equations are different, an additional transmission 
condition is needed. This additional condition can be chosen in two different ways. At the boundary of 
the micropolar material one can prescribe the value of the reduced couple traction , or, alternatively, 
the value of the rotation ipf (note that ip 2 = Lp 2 follows immediately from (|4*5jh ). We analyse these two 
cases separately in the next subsections. 

3.1 'Couple' transmission conditions 

The additional transmission condition in this case takes the form 

q+ = for 6 = 0. (49) 

Bearing in mind that, in this case, we need to match the solutions of governing equations having 
different orders, we use here a two-terms asymptotic analysis, so that a solution is searched for in the 
form 

w+(r,0) =r x F 1 (6)+r x+2 F 2 (0) + O(r x+4 ), r -> 0, (50) 
w-(r,6)=r x H 1 (6)+r x+2 H 2 (9)+0(r x+4 ), r -> 0. (51) 
Substituting ([50)) in ([4l))) i and collecting like powers of r, we obtain 

A 2 (r A Fi) = 0, A 2 (r A+2 F 2 ) = ^A(r A F : ). (52) 
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For A 7^ 0,1,2, the system ([5"2"j) admits the solution 

Ft = B^ sin(A6>) + B ( 2 1] cos(A6>) + sin[(A - 2)9) + b[ 1] cos[(A - 2)6), 

F 2 = B[ 2) sin[(A + 2)6»] + B (2) cos[(A + 2)9) + B {2) sin(A0) + B {2} cos(A6>) (53) 
B W o(l) 

+ 4AP Sin[(A " 2)6] + IaV C0S[(A _ 2)91 

Substituting ([5"Tj) in (|4^|) 9 and collecting like powers of r, we obtain 

A(r x Hi) = 0, A(r x+2 H 2 ) = 0. (54) 

For A ^ 0, 1,2, the system (|5"4")) admits the solution 

i?i = A^ ] sinXe + Af> cos A0, H 2 = A {2) sin[(A + 2)6) + A ( 2 2) cos[(A + 2)9). (55) 

Substituting the two-terms asymptotics with F1.2 and H\ t2 given by (p)3")) and ([53]) respectively, 
into the boundary and transmission conditions, we obtain a 12x12 homogeneous algebraic system, whose 
characteristic equation is 

(1 + rj) cos 2 (ttA) sin 4 (7rA) = 0, (56) 

so that the exponent A admits the values A = fc/2, where k = 1, 3, 5. 

Note that for A = 1/2 the first term in (|51[) for the classical elastic material corresponds to bounded 
elastic energy, while the first term in (|50|) for the micropolar material corresponds to infinite elastic 
energy. However, from the analysis of all boundary and transmission conditions, it is found that F\ (9) = 
in this case, so that the energetic requirements are fulfilled. 

For the special cases A = 0, 1, 2 the analysis is straightforward and the final asymptotic representation 
of the solution taking into account all terms A = fc/2, fc = 0, 1, 2, 3, is given by 

36* 56* 

w~ (r, 6) = ao + r 1 ^ 2 ai sin — h ra 2 cos 8 + r 3 ^ 2 a 3 sin h r 2 a 4 cos 29 + r b ^ 2 a% sin hO(r 3 ), r — > 0, 



w + (r,9) = a + r(a 2 cos6> + /?i sin0) + r 2 I a 4 ( — — + cos 29 ) + f3 2 sin 29 



2 2 

2G_ f 3-5/7 . 50 01 ^, 3 x 
01 3G + ^(3-^) { 5(lTij) Sm T + Sm 2 } + ° {r } ' r - °' 

The asymptotics of stress in the lower half-plane, occupied by the classic elastic material, as r — > are 
given by 

1 $ 3 (9 5 36* 

cr^ = --r~ 1/2 aiG- sin - + a 2 G- + -r 1/2 a 3 G- sin - + 2ra 4 G_ cos0 + -r 3/2 a 5 G_ sin — + 0(r 2 ), 

_ Z Z Zj Z Z 

a 23 = -r-V 2 axG- cos - + -r 1 / 2 a 3 G_ cos - - 2ra 4 G_ sin<9 + -r 3 / 2 a 5 G_ cos — + 0(r 2 ). (58) 
2 2 2 2 2 2 

The asymptotics of symmetric stress, couple-stress and skew-symmetric stress in the upper half-plane, 
occupied by the micropolar material, as r — > are given by 

a+3 = a 2 G + + 2rG + (a 4 cos9 + ft sinfl) - ? r 3 / 2 axG_i^^__M^ sin * + o^), 
a+ = PiG + + 2rG + (a 4VS m9 + (3 2 cc«9) - \r^a,G_ 1 ^ ' cos 5 + 0(r 2 ), 



(57) 
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M + = 2faG + l 2 (l + n) - ^VG„ + + cos i + 0(r), 

/4i = 7 «iG_ - — — ^ cos - + 0(r), (59) 

+ n ^ ,2/^ 2\ 1/? ^ 3 - 5?/ - (1 + 77) cos 6* . 9 

l4n = -2a 4 G+r 1 - rf) + r 1/2 ai G- '— i ^ sm - + Or, 

3 — 77 2 



^3 = -^f^ S ini + 



3-77 2 



The condition (|13|) becomes now /-i^ 2 ( 7 ' = 0, = 7r) = 0, which requires fa to vanish, except for the case 
where a tangential line load is applied along the crack edge. 

It is now possible to analyse the jump of the rotation component ip± across the interface and it is 
found that there is a mismatch between the micropolar material and the elastic one given by 

M = ~a 1 r- 1 ^ + ^ 1 - ^a^+far + O^' 2 ), r 0. (60) 

It is concluded that in the classical elastic material the solution shows a square-root singularity, as 
in the classic case. Moreover, in the micropolar material, the symmetric stress and the couple-stress are 
bounded, and only the skew-symmetric stress shows a square-root singularity. 

The energy release rate, computed in the same manner as in Sec. 12.31 is given by 

J* = l G - a l ( 61 ) 
which shows that the constant a± plays the role of stress intensity factor. 

3.2 'Rotation' transmission conditions 



The problem is defined by the equations (|46j) with the boundary and transmission conditions (|47[) and 
l|48p. The additional transmission condition used here is the continuity of the rotation vector, namely 

<pt = tpY for 9 = 0. (62) 

We search for a solution again in the form ([50)) and (|5"Tj) . For A 7^ 0, 1, 2, the characteristic equation now 
takes the form 

(1 + ?/) 2 sin 2 (7rA)[l + 77- (77- 3)cos2ttA] 2 = 0, (63) 
so that the exponent A admits real values, namely integer positive numbers and 

1 77 + 1 

A = ± — arccos h k, (64) 

2n rj — 3 

where k is a non negative integer. 

In the particular case r\ = 1, the complete asymptotics of the solution up to the forth order is given 

by 

w (r, 9) = a Q + ra 1 cos 9 + r 2 a 2 cos 29 + r 3 a 3 eos30 + r 4 a 4 cos 49 + 0(r 5 ), r -> 0, (65) 
w + (r, 9) — a + rai cos9 + r 2 a 2 + r 5 ^ 2 fa (cos ^— — cos - ] + r 3 a 3 ( ^ cos 6* — ^ cos 39 



- r 7/2/ 



fa Uos y - cos y ) + r 4 {(^p - 04) cos 49 + (2a 4 - g|) cos 20 + ^} + 0(r 9 / 2 ), r -> 0. 
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The asymptotics of stress in the classic elastic material occupying the lower half-plane as r — > 
become 

<rf 3 = aiG_ + 2ra 2 G- cos0 + 3r 2 a 3 G_ cos26> + 4r 3 a 4 G_ cos36> + 0(r 4 ), 

(66) 

a 23 = -2ra 2 G- sin (9 - 3r 2 a 3 G_ sin 29 - 4r 3 a 4 G_ sin 36* + 0(r 4 ), 

Finally, the asymptotics of symmetric stress, couple-stress and skew-symmetric stress in the microp- 
olar material occupying the upper half-plane as r — > become 

a+ 3 = o>iG + + 2m 2 G+ cos6» - 3r 3/2 /?iG + sin 9 sin - + 3r 2 a 3 G+ 

f 9 39\ 

-10r 5/2 /? 2 G+sin 2 - (2cos-+cos y J + 4r 3 a 4 G+ cos<9(2 - 2 cos (9) + 0(r 7/2 ), 

1 / 36\ 

a+ 3 = 2ra 2 G+sm9- -r 3/2 fi 1 G + f 3 sin - + 7sin — j + 3r 2 a 3 G+ sin 29 

A r s/*P 2 G + ( 5 sill 6 _ 9 sin ^ + J-r 3 G+ sin 0[a 2 - (a 2 - 18a 4 ^ 2 ) cos 20] + 0(r 7 / 2 ), 

Z \ Z 2 J 6t z 



w 



(67) 



o / q oQ \ 5 / Q 

At^i = --r 1/2 /?iG+Z 2 5 sin - + sin— + 12ra 3 G+Z 2 sin0 - -r 3/2 /3 2 G+Z 2 3 sin - + 7 sin 
+24r 2 a 4 G + Z 2 sin26' + 0(r 5/2 ), 

& = AI+ = - ^r^G+Z 2 ( 5 cos J _ cos |) + \r 3 / 2 h G + l 2 ( 3 cos 2 - 7 cos |) 
+4r 2 G+ sin 2 6>(a 2 - 12a 4 / 2 ) + 0(r 5 / 2 ), 

r+ = -r" 1/2 /3i G + Z 2 cos - - 6a 3 G+Z 2 + y r 1/2 /3 2 G+Z 2 cos - - 24ra 4 G+Z 2 cos 9 + 0(r 3/2 ), 

r+ = ^r- 1/2 /3iG+Z 2 sin^ - y r 1/2 /3 2 G+Z 2 sin - - 4rG+(a 2 - 6a 4 Z 2 ) sin 9 + 0(r 3/2 ). 

In this case, the reduced couple traction q± along the interface has the form: 

qf{r, 9 = 0) = -Q^^^G+l 2 - Vdr 3 ' 2 p 2 G + l 2 + 0(r 5 / 2 ), r -> 0. (68) 

Of course, it is possible to construct the asymptotic solution for arbitrary \r)\ < 1. However, the 
general form for arbitrary r\ is rather complicated. As an additional example, we provide the result for 
77 = 1/3 and up to the first four terms (A = 0, 1/3, 2/3, 1): 



w (r, 9) = a$ + ra,\ cos 9 + r 2 a 2 cos 29 + r 7 ^ 3 a^ ^ 
+r 8/3 a 4 (sin^ ]= cos ^) + 0(r 3 ), r^O 



76> 1 79 

a 3 I sin = cos — 

3 V3 3 



V3 

w + ( r , 9) = a + rai cos + r 2 a 2 Q + i cos 26>^ ( 69 ) 



7/3 / 10 . 76* 10 19 1 . 9 7 

+r ' as — sin = cos sin — I = cos - 

V 9 3 3^3 39 3 3^3 3 

8/3 /ll .86> 11 86» 8 . 26> 8 26>\ n . 
+r/« 4 ^ y sm y + y 7! cos y --sm y -— cos y j+G(r ), r ^ 0. 

It is found that the behaviour of the displacement and stress fields is similar for any — 1 < r\ < 1. 
In particular, stresses are always bounded in the classical elastic material, while in the micropolar 
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material singular behaviour appears only in the skew-symmetric stress (with different level of singularity 
depending on the value of if). Moreover, the energy release rate is always zero for any — 1 < r] < 1. This 
shows that this type of transmission conditions does not allow for the propagation of the crack along 
the interface, and thus has limited physical meaning. 



3.3 The special case r/_ = — 1 and 77+7^—1 

In the particular case of r/_ = — 1 and rj+ =/= —1, the governing equation fjl 1[1 and the traction- free crack 
face conditions (|T2"j) for the material in the lower half-plane are satisfied by the classical solution defined 
by the field equation 

Aw~ = 0, (70) 



together with the following boundary condition 



dx 2 



for xi < 0, x 2 = 0. (71) 



Moreover, in this case the couple-stress and skew-symmetric stress fields in the lower half-plane identically 
vanish, so that the transmission conditions along the interface x\ > 0, x 2 = become 

w + =w~, pj = 023, ( 72 ) 
together with one of the following additional conditions 

qf=0, iff = ^, (73) 
which correspond to the cases investigated in Sees. 13.11 and 13.21 respectively. 



4 Discussion and conclusions 

In the present work, the effects of strain rotation gradients on a stationary Mode III crack along the 
interface between dissimilar couple-stress elastic materials have been analytically investigated by per- 
forming an asymptotic analysis of the crack tip fields. It is shown that solutions without oscillations 
appear in the following two cases: when the two materials are the same (homogeneous material) and 

when the two materials are dissimilar but 77^ = 77 = 1 . In the latter case, the solution displays the 

same r~ 3 / 2 singularity (appearing in the skew-symmetric stress components) as for the problem of a 
crack in an homogeneous CS material. In other cases, the solution exhibits oscillatory behaviour in the 
vicinity of the crack tip, with the overlapping zone becoming increasingly large as the ratio rj between 
the characteristic lengths in one of the materials approaches the value —1. The energy release rate has 
been calculated by means of the conservation J-integral. It is shown that contributions of the integrals 
along the crack faces have to be retained and the additional boundary condition along the crack edge 
(Koiter, 1964) is essential to guarantee that the generalized J-integral (Lubarda and Markenscoff, 2000) 
remains bounded. This additional boundary condition has always been omitted in the earlier literature 
because so far only the symmetrical problem in homogeneous materials was discussed. The boundary 
condition along the crack edge breaks the symmetry and becomes fundamental for interface problems. 

The special problem of a crack along the interface between couple-stress and classical elastic materials 
has also been addressed. Two types of transmission conditions have been considered: 'couple' and 
'rotation' transmission conditions. In the former case, it is assumed that the couple-stress traction is 
continuous, and thus vanishes, at the interface. In the latter, it is assumed instead that the rotations 
are continuous at the interface. It turns out that the solutions are quite different in the two cases and 
it is not possible to satisfy simultaneously both type of transmission conditions, so that a mismatch is 
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always present at the interface, resulting either in a non-balanced couple-stress or a discontinuity in the 
micro- rotations. It is shown also that the special case rj_ = — 1 and n + ^ — 1 reduces to the problem of 
a crack at the interface between classical and couple-stress elastic materials. 
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